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ABSTRACT 

We present a method for computing the net transmission of X-rays emitted by shock-heated plasma 
distributed throughout a partially optically thick stellar wind from a massive star. We find the trans- 
mission by an exact integration of the formal solution, assuming the emitting plasma and absorbing 
plasma are mixed at a constant mass ratio above some minimum radius, below which there is as- 
sumed to be no emission. This model is more realistic than either the slab absorption associated with 
a corona at the base of the wind or the exospheric approximation that assumes all observed X-rays are 
emitted without attenuation from above the radius of optical depth unity. Our model is implemented 
in XSPEC as a pre-calculated table that can be coupled to a user-defined table of the wavelength 
dependent wind opacity. We provide a default wind opacity model that is more representative of real 
wind opacities than the commonly used neutral ISM tabulation. Preliminary modeling of Chandra 
grating data indicates that the X-ray hardness trend of OB stars with spectral subtype can largely be . 
understood as a wind absorption effect. 

Subject headings: stars: early type — stars: winds, outflows — stars: mass-loss — radiative transfer 


1. INTRODUCTION 

The absorption of soft X-rays by the powerful, 
radiation-driven winds of OB stars has long been recog- 
nized as a significant effect both on the X-rays observed 
from these stars and on the physical conditions in their 
winds. The soft X-ray emission observed in OB stars by 
Einstein and ROSAT implied only modest wind attenua- 
tion of the X-rays, and thus ruled out significant coronal 
emission as a source of the ubiquitous X-ray emission 
seen in these massive stars (Cassinelli Sz Olson 1979; 
Cassineili & Swank 1983; Hillier et ah 1993; MacFar- 
lane et al. 1994). For this and other reasons, the wind- 
shock paradigm for the production of X-rays in OB stars 
has become generally accepted (e.g., Owocki et al. 1988; 
Pauldrach et al. 1994; Feldmeler et ah 1997a,b; Kahn 
et al. 2001), although many aspects are still poorly un- 
derstood. 

Using X-ray observations to constrain X-ray produc- 
tion mechanisms requires proper account of the ab- 
sorption of distributed sources of X-ray emission. The 
amount and wavelength dependence of the wind ab- 
sorption can be used as a diagnostic of the loca- 
tion/distribution of the shock-heated plasma and of the 
wind mass-loss rate, especially in terms of its effect on in- 
dividual line profile shapes (Owocki & Cohen 2001, here- 
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after OCOl). Even simply deriving an intrinsic X-ray 
luminosity for energy budget considerations requires cor- 
rectly modeling the significant attenuation of the emit- 
ting X-rays, especially in the dense winds of O super- 
giants (Hillier et al. 1993: Owocki k. Cohen 1999). 

Because the emitting plasma is spatially distributed 
throughout the absorbing wind, simple prescriptions for 
the attenuation can be inaccurate. Specifically, the com- 
monly used slab model of absorption, appropriate for an 
intervening ISM cloud where all of the emission origi- 
nates beyond the absorbing medium, has transmission 
with an exponential dependence on slab optical depth, 
and strongly overestimates the amount of attenuation us 
the wind becomes optically thick. Furthermore, bemuse 
hydrogen and helium are ionized, models for neutral gas 
significantly overestimate wind opacity. Unfortunately, 
due to the lack of easily available tools for incorporating 
appropriate radiative transfer and opacities, inadequate 
models, such. as those intended for neutral interstellar 
medium absorption, are routinely applied to account for 
wind attenuation. This is in spite of the fact that a num- 
ber of previous works have recognized the necessity for 
and applied appropriate wind absorption prescriptions 
(Hillier et al. 1993; Pauldrach et al. 2001: Owocki <k* Co- 
hen 2001: Oskinova et al. 2006). 

We have developed a method for implementing an ex- 
act solution to a realistic model of the radiation transport 
that can be easily combined with a pre-calculated opacity 
table to find the wavelength dependent emergent. X-ray 
flux from a stellar wind. We also provide a reasonable 
default opacity model that can be used for most OB star 
winds. Our analysis tool, which wo name wivdiab^. for 
mind ^ able absorption, can be combined with an indepen- 
dent plasma emission model, such as the As trophy si cal 
Plasma Emission Code (apec) (Smith et al. 2001) that is 
widely used in fitting stellar X-ray spectra. This can be 
used to realistically model the low-resolution CCD spec- 
tra that are produced in large quantities by surveys of 
clusters and OB associations with Chandra and XMM- 
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Newton (e.g., Sana et al. 2006; Wang et al. 2008). it 
can also be used to model grating spectra in detail, and 
provides a means of disentangling the wind absorption 
effects from the emission temperature effects that both 
appear to contribute to the recently discovered trend in 
the morphology of OB star spectra observed at high res- 
olution with the Chandra gratings (Walborn et al. 2009). 


2. RADIATION TRANSPORT MODEL 

In this section, we derive an expression for the frac- 
tion of X-rays transmitted from a massive star wind 
as a function of the wavelength-dependent opacity, the 
mass-loss rate, and the wind velocity law. We make 
assumptions similar to those made in Owocki &: Co- 
hen (2001): we model the wind as a spherically sym- 
metric two-component medium, where a small fraction 
of the wind is heated to X-ray emitting temperatures 
(Tx ~ 1 — 10 MK), while the bulk of the wind is com- 
posed of relatively cool material (T w \ n d/T e ff ~ 0.5 — 1) 
that can absorb the X-rays via the bound-free opacity of 
the moderately ionized metals. We assume that the X- 
ray emission turns on at some radius Ro > E*. We also 
assume that the temperature distribution of the X-ray 
emitting plasma is the same over the the entire emitting 
volume. 

The assumptions we make regarding the distribution 
of the X-ray emitting plasma are footed in the available 
observational evidence, as well as in the results of ex- 
tensive theoretical simulations. Detailed studies of emis- 
sion line profiles as well as constraints from forbidden- 
to-inter combination line ratios in He-like ions support a 
picture in which X-ray emission starts at Ro ~ 1.5 R* for 
all .observable ions (Cohen et al. 2006; Leutenegger et al. 
2006; Cohen et al. 2010), with a roughly constant fill- 
ing factor above the onset radius. Numerical simulations 
have typically found onset radii for strong shocks that 
are comparable to this (Owocki et al. 1988). Feldmeier 
et al. (1997b) performed simulations of winds seeded with 
base perturbations. They have found that clump-clump 
collisions are important, and that the resulting X-ray 
emission is distributed primarily below ten stellar radii. 
Runacres & Owocki (2002) have also found X-ray emis- 
sion with onset radii of order 1.5 stellar radii and extend- 
ing out to large radii. However, the results for the onset 
radius are also sensitive to the treatment of the scattered 
radiation field (Owocki & Puls 1996; Owocki 2009). 

The observed X-ray luminosity as a function of wave- 
length is given by 

L x = Atc J dV Vx (r)e-^- x \ ( 1 ) 

where rj\(r) is the X-ray emissivity, and r(r,^, A) is the 
continuum optical depth of the dominant cool component 
along a ray from the emitting volume element to the 
observer. 

The optical depth can be derived as in Owocki & Cohen 
(2001) for a smooth wind. It is given by the integral 



dz* /c(A) p(r x ) . 


( 2 ) 


Here p and z are ray coordinates, with impact parameter 
p — \J 1 - p 2 r and distance along the ray : = pr. where 
f.t is the direction cosine to the observer at local radius 


r. k(A) is the atomic opacity of the wind, and /;(/•) is 
the density of the wind. Using the continuity equation. 
p(r) = A//47rr 2 u(r), where M is the stellar mass loss 
rate, and defining the characteristic wind optical depth. 


we can write 
where 
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is a dimensionless integral that depends purely on (ho 
ray geometry. Here R . * is the stellar radius, is the 
wind terminal velocity, and w(r) = v{r)/v^ is the sealed 
wind velocity. Note that we have assumed that a(A) is 
constant throughout the wind; we further discuss this 
assumption in § 3. We take the velocity to follow a beta 
law: v = t>oc(l - R*/r) 0 . We also take 9 — 1 in this 
paper as a good approximation for many O star winds: 
however, evaluation for general values of 3 is not difficult. 

The emissivity is assumed to scale with density 
squared, as in Owocki & Cohen (2001). We ignore the 
Doppler shift of the emitted X-rays and write 


r)\(r > Ro) = 7)\$ ? ^ . (6) 

Po 

This expression gives the correct radial dependence for 
the total line strength, but discards information about 
the profile shape. This is a justified approximation in cal- 
culating the broadband X-ray transmission of the wind. 
Here we assume that X-ray emission begins at a mini- 
mum radius R 0 , with po = p(R o) and t/a.o = *m(/?u). In 
this paper, we will assume that the X-ray filling factor 
is constant with radius; it is trivial to add a power-law 
radial dependence, as in OC01. We also assume that X- 
ray emissivity follows the same radial distribution at all 
observable wavelengths. 

The model described in the preceding paragraphs is 
illustrated in Figures 1 and 2. The model emissiviu and 
transmission are visualized separately in Figure 1. and 
together in the left panel of Figure 2. The right pane! 
of Figure 2 gives the net transmission for the exospheric 
approximation for comparison, w'hich we discuss at more 
length below. 

Using Eqs. 1-6, we can calculate the transmission of 
the wind as a function of r*. The transmission is simply 
the observed flux (Eq. 1) divided by the unattenuated 
flux, which can be found by setting r m = 0 in the same 
equation: 


_ L X ( T.) _ fdVp*e- 
L x ( 0) fdVp* • 


To account for the attenuation in the numerator, it is 
convenient to define an angle- averaged transmission from 
each radius r. 
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Figure 1. This diagram shows the X-ray emission and absorption 
properties of a model stellar wind. The observer is on the left. The 
black disc denotes the star, and the thin black line shows the onset 
radius for X-ray emission (1.5 Ft • in this case). Contours of con- 
stant X-ray emissivity (proportional to density squared) are shown 
with dashed lines at intervals of an order of magnitude in differen- 
tial emissivity. The absolute scale is normalized to the maximum 
emissivity at the onset radius of X-ray emission. Contours of con- 
stant continuum optical depth calculated for r, = 3 are shown with 
solid lines. 



Figure 2. The left panel is the same as Figure 1, but plotting the 
net emissivity (product of intrinsic emissivity times transmission, 
including occultation by the stellar core) instead of the intrinsic 
emissivity; this net emissivity has been renormalized so that the 
point of maximum net emissivity has a value of unity. The right 
panel shows the same plot in the exospheric approximation, where 
the grey zone is inside the radius of (radial) optical depth unity 
and thus does not contribute to the observed X-ray flux. 

where 



gives the p coordinate of occultation by the the stellar 
core. Some X-rays are obscured by the stellar core even 
when the wind is transparent: 

r 0 (r)-?M)=i(l+M.). (10) 

Integrating over shells at all radii, the net transmission 


is thus 

r(r*) = 

We can further evaluate this expression by substituting 
the continuity equation, and by defining the inverse ra- 
dial coordinate u = R m /r : 


Jn 0 cfrrV(r)r(r,r.) 

IZ <*' r V( f ) 
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T(r.) 


' Jo'' dvw~ 2 (u)T{u. r») 
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( 12 ) 


Figure 3 shows the angle-averaged transmission T(u) 
for selected values of r *. Figure 4 shows the following 
quantity: 


Lx(u,t.) fidu'w-*(u')T(u',r.) • 

Lx(u 0 ,t*) Jq° du , w~ 2 (u')T(u',Ti,) 

This can be thought of as the cumulative distribution 
of observed X-ray emission, starting from u = 0 (ft = 
oo) and integrating in to u = u 0 . Together, figures 3 
and 4 show the relative importance of transmission and 
emission as a function of radius in stellar winds. For the 
entire range of interest in r*, essentially all of the wind 
down to uo contributes to the observed X-ray flux. 

Figure 5 shows T(r*) for this model, along with emu- 
parisons to two other absorption prescriptions: a simple 
intervening absorber. T = e~ r . appropriate for a coro- 
nal slab model, e.g. as implemented in the XHPEC mod- 
els wabs or tbabs: and an exospheric approximation (e.g. 
Owocki & Cohen 1999), where T = 0 below the radius 
of optical depth unity, and T = 1 everywhere above it: 


T e xo(7* ) 


/ 0 U " duu'--{u) 
Jo” dvu- 2 (u) ‘ 


(H) 


where u x — min(ui(r„), uq). The inverse radial coordi- 
nate of optical depth unity is given by evaluating the 
optical depth integral (Eq. 5) along a radial ray (p = 0, 
z = r), with the result 


= 1-e i/Tm , (15) 


for 0 = 1. Note that in Figures 3, 4 and 5 we have used 
Q ~ 1 and.fto ~ 1.5ft*. We stress that Eq. 14 is simplv 
the consequence of using a step function for ihe artglc- 
averaged transmission T in Eq. 12. rather than the more 
realistic expression given in Eq. 8. 

The transmission of our model falls off much more 
gradually than e“ r (Figure 5), but it is also more ac- 
curate than the exospheric approximation, especially at 
moderate optical depth. The exospheric approximation 
has the correct asymptotic behavior for large values ofr*. 
but overestimates the transmission by a fixed factor. Tlu* 
exospheric transmission at large r, can be brought into 
agreement with wmdtabs by multiplying the exospheric 
r m by three. 

It would be possible to generalize the radiation trans- 
port model described in this section to include poros- 
ity by introducing an appropriate definition of effective 
opacity, as in Oskinova et al. (2006) or Owocki & Cohen 
(2006). However, a specific implementation of this and 
discussion of its consequences is left to a future study. 
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Figure 3. Angle averaged transmission T as a function of shell 
inverse radial coordinate u — R m / r, with different curves for dif- 
ferent characteristic wind optical depths r*. The equivalent plot 
in the exospheric approximation would be a step function at uj, 
the inverse radial coordinate of optical depth unity. Note that the 
transmission is less than unity even for r- = 0 due to occultation 
by the stellar core. 



Figure 4. Fraction of X-ray emission originating from outside 
inverse radius u, normalized to total X-ray emission integrated to 
uq = 2/3. For all but the very optically thick case of r- = 10. the 
emission comes from a wide range in u . 


3. OPACITY MODEL 

To determine the wind transmission as a function of 
wavelength, we need to account for the wavelength de- 
pendence of the optical depth through the opacity, which 
we write as 

r*(\) = k(A)E* , (16) 

where 




M 

47Tft„1/oo 


(17) 


is the characteristic mass column density of the wind (in 
g cm -2 ). 

The continuum opacity of a stellar wind in the X-ray 
band can be calculated by summing the contributions of 
each constituent species. Thus, we must know the atomic 


Figure 5. Comparison of transmission of three different mod- 
els: coronal slab (e”" T ), exospheric, and more realistic wind model 
(windiabs). The fixed parameters are 0 = 3. Rq = 1.5. 

cross-sections, ionization fractions, and elemental abun- 
dances. Of these three, the cross sections are known to 
sufficient accuracy (e.g. Verner k Yakovlev 1995): the 
ionization balance may contribute some uncertainty in 
the calculation of the opacity, but is usually not the dom- 
inant source of error; and uncertainties in the elemental 
abundances are typically the most important. 

The opacity due to photoionization of a given shell of 
any individual species scales approximately as K t (A) x A ;i 
above the threshold energy £ t]) of the shell. Because mul- 
tiple species are usually important for the X-ray opac- 
ity of astrophysical gas, the run of opacity with wave- 
length has a characteristic sawtooth shape, with individ- 
ual teeth at the ionization threshold energies of dominant 
ionization stages of abundant elements. 

O star winds are photoionized. with H and He fully 
stripped (although He may recombine in some dense 
winds; see below), and most other elements mainlv in 
charge states 4-3 and 4-4. Thus, the opacity of stellar 
winds in the range 1 A < A < 40 A is dominated by 
K-shell absorption in C, N. and O, since they are the 
most abundant elements. Significant contributions from 
K-shell absorption in Ne, Mg. and Si are also present, hs 
well as Fe L-shell absorption. 

The opacity of adjacent ionization stages of the same 
element are usually comparable, with the exception that 
the photoionization threshold energy is shifted. The* ef- 
fect of a moderate shift in ionization balance on the opac- 
ity is relatively minor, since none of the strong X-ray 
emission lines observed from O star winds falls between 
the threshold energies of adjacent stages of abundant 
ions. However, the difference in threshold energies, and 
hence the broadband absorption, between neutral and 
wind material is significant. 

Therefore, while it is important to use an appropriate ' 
model for the wind ionization, it is sufficient for many ap- 
plications to use a single approximate ionization balance 
to model all O star winds. This is true even though (he 
ionization balance can vary to some extent with radius, 
and also is different in different stars. 

If it is desired to model the opacity of a particular star, 
it is possible to construct a detailed model opacitv for a 
stellar wind by using the output of a radiative transfer 
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model such as CMFGEN (Hillier et al. 1993; Hillier & 
Miller 1998; Cohen et al. 2010, Zsargo et al., in prepara- 
tion). 

To illustrate the importance of various assumptions in 
opacity modeling, Figures 6 and 7 compare several model 
opacities. Figure 6 compares neutral interstellar medium 
opacity and a simple model for the opacity of a stellar 
wind. Both assume solar abundances (Asplund et al. 
20.09). In the wind model we assume an ionization bal- 
ance with hydrogen and helium fully stripped, oxygen 
and nitrogen in the +3 charge state, and all other ele- 
ments in the +4 charge state. This is a good approxima- 
tion to the ionization balance in a typical O star wind, 
which is set by photoionization from the photospheric 
UV field. The model wind opacity is much lower than 
the model ISM opacity, especially at long wavelengths, 
mainly due to the ionization of hydrogen and helium. 
The shift in ionization threshold energies due to the pres- 
ence of more highly ionized species is also clear. 

Figure 7 shows three stellar wind model opacities, and 
thereby illustrates the relative importance of two effects 
on the opacity: the elemental abundances, and the ion- 
ization balance. The solid line gives the same solar abun- 
dance wind model described in the previous paragraph, 
while the dashed and dotted lines give models particu- 
lar to C Pup, using non-solar abundances specific to the 
star derived from detailed CMFGEN modeling (Bouret 
et ah, in preparation). The dashed line gives the ac- 
tual opacity at ~ 2 R+ from this CMFGEN model, while 
the dotted line uses the simplified ionization structure of 
the solar abundance wind opacity model. The fact that 
the realistic CMFGEN £ Pup model and the simplified 
version are so similar indicates that a simple ionization 
balance is typically adequate to describe wind opacity in 
many - cases, as long as the dominant ionization stages 
are relatively accurate. On the other hand, the differ- 
ence between the £ Pup models and the solar abundance 
model shows that the opacity model depends strongly on 
the abundances of the most common elements other than 
H and He (typically C, N, and 0). Note that the Bouret 
et al. abundances for C Pup are both non-solar in the ra- 
tio of CNO and also have sub-solar metallicity; it is the 
sub-solar metallicity that accounts for the lower opacity 
of the C Pup models at short wavelengths compared to 
the solar abundance model. 

We have made one important simplification in our 
modeling: in Eq. 4, and throughout this section, we have 
assumed that the opacity is independent of radius. As 
shown in Figure 7, moderate changes in wind ionization 
do not strongly affect the opacity, so in most cases this 
is a justified approximation. The important exception 
is the ionization of helium; in sufficiently dense winds, 
lle ++ may recombine to He + in the outer part of the 
wind, which greatly increases the opacity, especially at 
long wavelengths (Pauldrach 1987: Hillier et al. 1993). 
As long as the change in ionization occurs sufficiently far 
out in the wind, geometrical effects as described in § 2 
are not important, and the absorption due to He + can be 
treated as an additional slab between the X-ray emitting 
regions and the observer, i.e. using T = e~ r . 

4. MODEL IMPLEMENTATION 

The numerical evaluation of Eq. 12 is not prohibitively 
expensive, but it is typically not fast enough to allow its 



Figure 6. Comparison of a neutral interstellar absorption model 
(including dust grains), tbabs (Wilms et al. 2000); and an O star 
wind model with an assumed simple ionization structure. The 
opacity k(X) is scaled with A 3 to allow a better comparison of the 
two models at short wavelengths. Both models use solar abun- 
dances (Asplund et al. 2009). Note that even below the 0 3+ edge 
near 20 A, the wind opacity model is still about 40% lower than 
the ISM model, mainly due to the ionization of H and He in tlx* 
wind model. 



Figure 7. Comparison of three different model O star wind opaci- 
ties: the solar abundance O star wind model (with assumed simple 
ionization structure) shown in Fig 6 (solid); a CMFGEN model ap- 
propriate for £ Pup (dashed) using the abundances of Bouiet et al 
(in preparation): and a simplified version of this CMFGEN model 
(dotted: red in the elect ronic version) where the ionization balance 
is the same as in t he solar abundance model. Note tlmt the Hornet 
et al. abundances for £ Pup are subsolar as well as having alirrrd 
CNO abundance ratios. 

use in an automated spectral fitting routine, such as that 
in XSPEC12 (Arnaud 1996) or ISIS (Houck L Denicola 
2000). It is thus preferable to compute the transmission 
on a grid in r* for a given set of wind parameters. 

Given a tabulation of the model wind opacity, as de- 
scribed in § 3. in addition to the tabulation of 7'{rJ. 
one may then calculate the transmission as a function oi 
wavelength, T( A), with only one free parameter . the char- 
acteristic wind mass column density This parameter 
is analogous to the neutral hydrogen column density in a 
slab absorption model such as the XSPEC wabs or tbabs 
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Figure 8. Transmission as a function of wavelength for ionized 
wind absorption model ( windtabs> black) and for neutral slab ab- 
sorption (tbabs, grey; red in the electronic version). Three values of 
absorbing column are given; for windtabs, the degree of absorption 
is specified by the characteristic mass column density while for 
tbabs it is given simply by the mass column £. 

ISM absorption models (which are also sometimes used 
to approximate wind absorption). 

We have implemented this as a local model for XSPEC 
12 9 . A user can calculate T(r *) for a given set of pa- 
rameters (i.e. 0 t Rq) 9 with the results stored in a FITS 
! table. These implicit model parameters may be varied 
j by computing additional tables of T(r*) for each set of 
parameter values. However, the absorption model is not 
j very sensitive to these parameters over the range typi- 
cally inferred for winds of massive stars. The calcula- 
tion of T(r*) is controlled by a simple python script, and 
| computation of a table for a given set of parameters can 
be accomplished in several seconds on a modern work- 
station. The model opacity must also be supplied as a 
FITS file; different model opacities may be swapped in 
at run time. When windtabs is used in a spectral fitting 
program, model transmission is calculated as a function 
of energy or wavelength using the supplied FITS tables 
and the one free model parameter, E*. 

The elemental abundances, which enter into the trans- 
mission through their effect on the opacity, cannot be 
varied as fit parameters in our model. This is a choice 
we have made in the model implementation, both for 
computational ease and simplicity of user interface, and 
! because there is not enough information in X-ray spectra 
I to constrain elemental abundances through modeling of 
absorption of X-ray emission alone. Abundances should 
be inferred by other means, e.g. from global fits to UV 
I and optical spectra, or from fits to X-ray emission line 
l strengths . These abundances can be used to compute a 
new opacity table for a given star. 

Figure 8 gives the model transmission for windtabs us- 
ing three different values of XI*, and using our standard 
0 star wind opacity model. For comparison, this figure 
i also shows the transmission for the neutral absorption 
model tbabs for comparable mass column densities E. 

9 Source code is freely available under the General Pub- 
lic License and may be obtained by request from the cor- 
responding author, or from the subversion repository at 
https://windprofile.svn.sourceforge.net/svnroot/windprofile/trunk 


Note that E refers simply to a slab mass column density, 
while E* refers to a characteristic mass column density 
in the context of a stellar wind (see Eq. 17). 

5. DISCUSSION 

5.1. Advantages of windtabs over exospheric or 
erponential attenuation 

The exact method we have presented here for modeling 
the X-ray radiative transfer for a distributed source em- 
bedded within a stellar wind is crucial for analyzing the 
X-ray emission observed from O stars. Compared to the 
exponential, neutral slab absorption (excess over ISM) 
model that is usually employed, the windtabs model ac- 
curately reproduces the much more gradual decrease in 
transmission with increasing wind column density and 
opacity. This can be seen in Figure 8, where the exponen- 
tial transmission model shows an unrealistically sudden 
decrease in transmission when the wind becomes opti- 
cally thick. 

Because the opacity of the bulk wind is a relatively 
strong function of wavelength, the inaccuracy of the ex- 
ponential transmission model will lead to errors in the 
broadband spectral energy distribution of a model ap- 
plied to individual 0 stars, leading to misinterpretations 
of the associated emission model components. This ap- 
pears to be the case in the study of Chandra grating 
spectra in Zhekov & Palla (2007). which invokes excess 
exponential, neutral ISM absorption to account for the 
assumed wind attenuation. Most likely because the ex- 
ponential transmission radiation transport model, and 
also the neutral ISM opacity model, significantly overes- 
timates the degree of attenuation, these authors find ab- 
sorption beyond the ISM column for only one star, even 
though some additional wind attenuation is expected for 
most of the stars in their sample. Additionally, although 
the authors do not comment on it, their determination of 
elemental abundances for each of the stars shows a consis- 
tent trend of abundance correlated with the wavelength 
of the dominant line or lines from each element. Such 
an effect would be expected if the wavelength-dependent 
wind attenuation is not accurately accounted for. 

It has long been noted that the exponential attenua- 
tion treatment is not well suited to modeling OB star 
X-rays. Hillier et al. (1993) used an exact treatment in 
their modeling of C Pup, discarding the possibility of a 
coronal model on the basis of the strong soft observed 
X-ray flux of ( Pup, together with its relatively high 
mass- loss rate. Cohen et al. (1996) showed that an exo- 
spheric treatment, rather than an exponential treatment, 
is important for understanding the observed EUV and 
soft X-ray emission from the early B giant, e CMa. The 
exospheric approximation was used by Owocki k Cohen 
(1999) to model the effect of wind attenuation of X-rays 
in order to explain the observed L x /Lb 0 i ~ 10“ 7 rela- 
tionship and its breakdown in the early B spectral range 
where hot star winds become optically thin to X-rays. 
An exospheric treatment was also used by Oskinova et al. 
(2001) to analyze the variability of X-rays from optically 
thick WR winds. And the exospheric framework lot ms 
the basis for the "optical depth unity'* relationship for X- 
rays in 0 stars, where the forbidclen-to-iutercoinbination 
line ratios of helium-like ions are claimed to implv for- 
mation radii that track the optical depth unity radius 
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as a function of X-ray wavelength (Waldron & Gassinolli 
2007). ■ 

However, the exospheric treatment underestimates the 
attenuation of the wind. If data need to be analyzed with 
a high degree of accuracy, then a more realistic treatment 
of X-ray radiation transport through the wind must be 
used; one that takes the inherently non-spherically sym- 
metric nature of the problem into account. This is es- 
pecially true when the location of the X-ray emission 
emerging from the wind is important, as with the inter- 
pretation of f /i ratios. As Figure 2 shows, the relative 
contribution of the emergent X-ray flux from different 
wind regions is incorrectly estimated in the exospheric 
approach, and as Figure 4 shows, the emergent flux has 
significant contributions from a wide range of radii. In 
fact, a number of previous works have implemented ac- 
curate X-ray radiative transfer prescriptions (e.g., Paul- 
drach et al. 2001; Owoc.ki & Cohen 2001; Oskinova et al. 
2006); however, none of these works provide an imple- 
mentation that Is available for use in other contexts. 

The windtabs model we have introduced here is not 
only more accurate in terms of the radiation transport 
than slab or exopheric treatments, but it has two addi- 
tional advantages that recommend its adoption for rou- 
tine X-ray data analysis and modeling of O stars. First, 
it is easy to use and has only a single free parameter, 
the characteristic mass column density E*, from which a 
mass-loss rate can be readily extracted. And second, it 
incorporates a default wind opacity model that is signifi- 
cantly different from, and much more accurate than, the 
neutral ISM opacity models that are usually used. Ad- 
ditionally, alternate user-calculated opacity models are 
easy to incorporate. 

5.2. Implications for trends in X-ray hardness with 
spectral type 

One application of windtabs to the interpretation of 
X-ray spectral data is for the analysis of the X-ray spec- 
tral hardness trend vs. optical spectral subtype recently 
noted . in Chandra grating spectra by Walborn et al. 
(2009). They report “the progressive weakening of the 
higher ionization relative to the lower ionization X-ray 
lines with advancing spectral type, and the similarly de- 
creasing intensity ratios of the H-like to He-like lines of 
the <v ions/’ The correlation of the overall X-ray spec- 
tral hardness with spectral subtype described by Wal- 
born et al. appear to be a direct result of wavelength 
dependent absorption effects. A detailed analysis is in 
preparation, but here we show a single suite of models 
in which the only variable parameter is the characteris- 
tic wind column density, £*. A single emission model, 
combined with windtabs attenuation, reproduces the ob- 
served broadband trend quite well. 

Figure 9 shows X-ray grating spectra for seven O giants 
and supergiants in the left-hand column, with the earliest 
spectral subtype (03.5) on the top, and the latest (B0) 
on the bottom, as in Walborn et al. (2009). The later 
spectral subtypes clearly have more soft X-ray emission, 
although the earlier subtypes still have non-negligible 
long-wavelength (A > 15 A) emission. In the middle 
column we show a four-temperature APEC (Smith et al. 
2001) thermal equilibrium emission model. We have cho- 
sen kT ” 0.1. 0.2, 0.4, and 0.6 keV: the first three com- 
ponents having equal emission measures and the hottest 


Table 1 

Adopted stellar parameters 


Star 

Spectral type* 

W H b 


yv. d ‘ 

HD 150136 

03.5 If* + 06V 

0.76 

0.073 

3.6 

C Pup 

04 I(n)f 

0.01 

0.160 

7.2 

f Per 

*07.5 IH(n) ((f)) 

0.115 

0.017 

0.76 

r CMa 

09 II 

0.056 

0.013 

0.58 

S Ori 

09.5 II 4- B0. 5 III 

0.015 

0.011 

0.49 

C Ori 

09.7 lb 

0.03 

0.019 ' 

' 0.85 

c Ori 

B0 la 

0.03 

0.020 

0.89 


a Spectral types adopted by Walborn et al. (2009). 
b Interstellar medium column density (10 22 cm" 2 ) 
c Characteristic wind mass column density {g cm* 2 ) 
d Characteristic wind number column density equivalent 
to E* (10 22 cm- 2 ) 

one having half the emission measure of the others. The 
same APEC model is used for all seven stars (with variable 
overall normalization) and is multiplied by a windtabs 
model and a tbabs model (for neutral ISM attenuation). 
The column density of the tbabs model is fixed at the in- 
terstellar value taken from FVuscione et al. (1994), with 
the exception of HD 150136, for which we inferred the 
ISM column density from E(B-V) (Mafz-Apellaniz et al. 
2004; Martins Sz Plez 2006: Vuong et al. 2003). The char- 
acteristic mass column density E„ in windtabs is fixed at 
a value computed from the “cooking formula* theoreti- 
cal mass- loss rate computed by Vink et al. (2001), using 
the measured terminal velocity of Haser (1995) and the 
modeled radii of Martins et al. (2005). The standard 
solar abundance wind opacity model (solid line in Fig- 
ure 7) was used in windtabs , and the APEC model abun- 
dances were set to solar. There are no free parameters 
in these models, and the temperature distribution has 
not even been significantly optimized to match the data. 
The adopted parameters are listed in Table 1 . 

As the middle column of Figure 9 shows, the- sim- 
ple, universal emission model reproduces the broadband 
trend very well. TYends in individual line ratios genet ally 
cannot be reproduced only by accounting for the vary- 
ing attenuation, as pointed out by Walborn et al. (2009), 
but note that the Ne IX (13.5 A) to Ne X (12.1 A) ra- 
tio does indeed vary due only to differential attenuation 
among the earliest spectral subtypes. The right-hand 
column in Figure 9 shows the same emission and ISM 
attenuation models as in the middle column, but with 
excess exponential (neutral ISM) attenuation accounting 
for the wind absorption, again according to the wind col- 
umn densities expected from the adopted mass-loss rates, 
radii, and terminal velocities. The exponential attenua- 
tion trend seen in the right-hand column is too strong for 
the earliest spectral subtypes and too weak for the latest 
ones, where the different ISM column densities actually 
dominate the trend. 

The contrast between the windtabs and exponential 
models is quite stark,, and indicates that the more re- 
alistic models should generally be used when analyzing 
X-ray spectra, both high-resolution and broadband. It is 
also impressive how much of the observed spectral hard- 
ness trend is explained by wind attenuation, in the con- 
text of a realistic model. Not only do quantitative analy- 
ses of the suggested line ratio trends have to he 'evaluated, 
but a global spectral modeling that allows for both omis- 
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Figure 9. Left column: a sequence of Chandra spectra of 0 giants and supergiants from Walborn et al. (2009): middle column* multi- 
temperature thermal emission model with windtabs wind absorption model; right column: same model as middle column, but with tbabi 
neutral slab absorption model. 


sion temperature variations and wind attenuation vari- 
ations should be undertaken in order to disentangle the 
relative contributions of trends in emission and absorp- 
tion to the overall, observed trend in the spectral energy 
distributions. 

6. CONCLUSIONS 

We have presented an exact solution to the radiation 
transport of X-rays through a spherically symmetric, 
partially optically thick O star wind, and shown that 
it differs significantly from the commonly used slab ab- 
sorption and exospheric models. Specifically, the trans- 
mission falls off much more gradually as a function of 
fiducial optical depth in the windtabs model as compared 
to the exponential model, leading to more accurate as- 
sessments of wind column densities and mass-loss rates 
from fitting X-ray spectra. As one example of the utility 
of windtabs. we have shown that when this more accu- 
rate model is employed, differential wind absorption can 
explain most of the observed trend in OB star X-ray 
spectral hardness with spectral subtype, and even may 
explain some of the line ratio trend. 

The windtabs model has been implemented as a custom 
model in XSPEC, and is as easy to use as the various ISM 
absorption models, having only one free parameter. In 
addition to the significantly improved accuracy of the ra- 
diation transport, windtabs has several other advantages. 


It incorporates a default opacity model much more ap- 
propriate to stellar winds than the neutral element opac- 
ity model used in ISM attenuation codes. Users can eas- 
ily substitute their own custom-computed opacity mod- 
els. And the fitted mass column density parameter for 
windtabs allows for the user to extract a mass-loss rate 
directly from their fitting of X-ray spectra of OB stars. 
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